Time Reversal Symmetry Breaking States Near Grain Boundaries 
Between d-wave Superconductors 
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In this paper we study the order parameter and density of states near a grain boundary between 
two d x 2_ y 2 superconductors. We examine broken time reversal symmetry near the interface. In 
particular we show that, under suitable circumstances, time reversal symmetry must be broken even 
when the order parameter is purely d x 2_ y 2 everywhere in space. 
PACS numbers: 74.50. +r, 74.72.-h 



The d x 2_ y 2 order parameter, appropriate to the hole- 
doped oxide superconductors, preserves time-reversal 
symmetry (TRS) in the bulk. At surfaces and interfaces 
it is now known that time-reversal symmetry may be 
broken. NIS-tunneling experiments by Covington et al. 
jlj indicate a time reversal symmetry breaking (TRSB) 
state locally at surfaces. Fractional fluxes at corners 
of interfaces in inclusion experiments by Kirtley et al. jj| 
strongly indicate that TRS may also be broken at grain 
boundaries. In this paper we discuss the origin of 

the TRSB-state and constrast TRSB at low and high 
transmission interfaces. 

At a surface or interface with low transmission, TRSB 
can occur in the presence of subdominant pairing inter- 
action in channels other than the dominant d x 2_ y 2. In 
this case order parameters corresponding to those chan- 
nels can appear near the interface K|. For example, if a 
subdominant pairing interaction is present in the s-wave 
channel, then, under suitable conditions, the order pa- 
rameter near the interface can have a d ± is symmetry. 
The order parameter thus breaks TRS locally. A pre- 
requisite for the TRSB-state is substantial pair-breaking 
at the interface, i.e. the misorientation of the surface 
normal to crystal a-axis should be close to j (mod 7r/2). 

The above is in constrast to the case where there is 
a reasonably high transmission probability of electrons 
across the interface. In this case the subdominant pair- 
ing interaction is not necessary for TRSB at the inter- 
face. [^| jic| , [ll|l TRSB occurs even when the order param- 
eter is purely d x 2_ y 2 everywhere near the interface. The 
origin of this TRSB-state is a proximity effect and it 
arises because the minimum energy state for the inter- 
face corresponds to a state with a finite phase difference, 
<-^X — Xr ~ Xl, across the junction. Here xl an d Xr are 
the phases of the order parameter on cither side far away 
from the interface. Ax is other than an integral mul- 
tiple of 7r for the TRSB-state. In this case states with 
minimum total interface free energy occur in pairs re- 
lated by time-reversal: if Ax corresponds to a state with 
minimum energy, there is also a non-equivalent but de- 
generate state with Ax' = — Ax- An interface at its min- 
imum energy configuration will have TRS spontaneously 
broken. In contrast to the one discussed in the last para- 



graph, this is the more likely route to TRSB when the 
transmission probability across the interface is moderate 
to high and when the misorientation between the two 
superconductors are close to ir/4 (mod tt/2). ]nl[L2] 

As is well-known, Josephson effects occur in the pres- 
ence of an interface with finite transmission. At a general 
phase difference between the two superconductors, a dis- 
sipationless current, J x , can flow through the interface. 
If the order parameter itself does not break TRS, then 
the current across the interface is always zero for Ax be- 
ing an integral multiple of %. However, under appropriate 
conditions there can be additional values of Ax where the 
net current across the interface vanishes. Previously one 
of us pj has explained how this can occur for a pinhole 
junction by considering the sum of contributions to the 
current from different parts of the Fermi surface. These 
new states with zero net current occur as a combined 
result of non-sinusoidal current-phase relationships and 
sign changes of the order parameter. The J x — states 
correspond to states where the junction energy is at a rel- 
ative extrema as a function of the phase difference Ax- 
In particular it can be shown that the new J x = states, 
if they exist, correspond to energy minima, fa] 

In this paper we study a planar interface with uni- 
form transmission. In general a current flows through 
the interface. The corresponding states possess finite 
flow energy densities even far away from the interface. 
Here we focus on the set of states with J x = for which 
this contribution is absent. In these states the gradient 
of the phase of the d-wave order parameter vanishes as 
x — > ±oo. We shall show that many of the statements 
concerning the new energy miminum states mentioned 
above for the pinhole || are still correct for the planar 
interface, provided appropriate minor modifications are 
made. If the order parameter is purely d x 2_ y 2, it is easy 
to verify that states which correspond to J x — with 
Ax — or 7T with x piecewise constant are always pos- 
sible. At not too small transmission across the interface, 
there may be other states with different Ax which also 
correspond to J x = and under appropriate conditions 
states with Ax ^ or tt will correspond to the mini- 
mum energy. We shall compare the free energies, order 
parameters and densities of states (DOS) of these J x = 
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states. Apart from its intrinsic interest, we shall see that 
the DOS provides an alternative view of the mechanisms 
for TRSB. A signature of a TRSB-state is that the ZEBS 
at s = are shifted away from the midgap and that spon- 
taneous currents along the interface are nucleated. MM 
The occurrence of zero energy bound states (ZEBS) 
for non-transmitting surfaces have already been exten- 
sively investigated || (and references therein) . For order 
parameters real up to a gauge transformation, ZEBS are 
present for the quasiparticle paths along which a sign 
change of the order parameter occurs. ZEBS are also 
common for interfaces with finite transmission if TRS 
is preserved (Ax = or tt). One can show rigorously 
fyl that ZEBS are present irrespective of the value of 
the transmission coefficient whenever there are quasipar- 
ticle paths such that a quasiparticle experiences a sign 
change of the order parameter both if it is transmit- 
ted or reflected. For interfaces between superconductors 
with large misorientation and in states which preserve 
TRS, ZEBS occur over a large part of the Fermi surface. 
The existence of these low energy bound states corre- 
sponds to severe pair-breaking near the interface. These 
ZEBS can be pushed to finite energies by allowing a fi- 
nite phase difference between the two superconductors. 
Correspondingly we shall show that the magnitude of the 
order parameter for the TRSB state (denoted simply by 
Ax 7^ below) is larger than the corresponding states 
with Ax = or tt. The formation of ZEBS and the sup- 
pression of the order parameter near the interface suggest 
that the Ax = or tt states are energetically unfavorable 
compared with the TRSB state. |l4| This is verified by a 
calculation of the free energy. 

TABLE I. Free energy per unit surface area calculated as 
for the junctions shown. The unit of 
. Nf is the normal state DOS. 
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For definiteness, we model the interface as an ideal, 
smooth barrier with a delta function potential. In this 
case, the interface can be parameterized by T> Ql the coef- 
ficient of transmission for normal incidence. The trans- 
mission coefficient T>((f>) for momenta pf at an angle <p 
with respect to the interface normal is given by 



V Q cos 2 



l - P shr 



(1) 



tions at the interface parameterized by 2?(</>) (see Refs. 

and jl2| for details). The numerical procedure is 
as follows: We start with an initial ansatz of the order 
parameter A(x), which in general possesses a phase dif- 
ference far away from the interface. For the given A(x), 
we obtain the quasiclassical green's function g(pf,e,x) 
along the imaginary axis (here e = ie n where e„'s are 
the Matasubara frequencies) via the method described 
in detail in No iteration is needed in this step. The 
correction to A(x) is determined by the weak-coupling 
gap equation involving the off-diagonal parts of g. At 
each iteration the current across the interface is also cal- 
culated. Then a gauge transformation depending on the 
calculated current is performed on the order parameter 
to relax A towards the state with J x = with the desired 
sign of J^x ■ This procedure is repeated until A(x) is con- 
sistent with the calculated g and the condition J x = 0. 
Note that once self-consistency is achieved, particle con- 
servation will be respected (see, e.g., Jll|). J x is then x 
independent and thus J x = at all x. After obtaining 
the self-consistent order parameter we evaluate the free 
energy via the " A trick" . |r| 

With the self consistent order parameter the green's 
function g was evaluated on the real energy axis again via 
the method in E3]. From this g we obtain the density 
of states. All the DOS below are obtained at energies 
e + i"f, with 7 = 0.05T C . This small imaginary part of 
the energy simulates a broadening of energy levels that 
would occur naturally in non-ideal systems. 

We first confine ourselves to pure d x 2_ y 2 order param- 
eter. We write 



A(p/,x) = 7] d {x)V2cos(2((j)- a)), 



(2) 



The quasiclassical green's function g obeys the usual 
equation involving A in the bulk, and boundary condi- 



which defines the complex order parameter r}d(x) with 
a = ax or ctR for the left and right side of the interface re- 
spectively. Here and an denote the orientations of the 
crystals on the two sides of the interface. They specify 
the angle between the a axis and hence the positive lobe 
of the order parameter with respect to the normal to the 
interface. We find that TRSB is most significant at low 
temperatures and when the misorientation 9 = olr — Ol 
is close to 7r/4 (mod it/2). 

As representative example we consider oil = and 
aji = tt/A at a relatively low temperature, T = 0.2T C . 
The order parameters are as shown in Fig [j] for both the 
states with Ax = and the ones corresponding to energy 
minima with Ax ^ 0. As can be seen from an examina- 
tion of Fig. [I] the phase difference Ax of the minimum 
energy state is ir/2. This state is degenerate with its time 
reversed partner —ir/2. The states with Ax = are also 
degenerate with the corresponding ones with Ax = tt 
with the same DOS. As claimed the order parameter of 
the Ax 7^ state has a larger amplitude than the one 
with zero phase difference for a given transparency. This 
difference decreases as V a decreases. At T> = 0.3 the dif- 
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ference in magnitudes is almost undetectable. The corre- 
sponding DOS are shown in Fig. |] for the two different 
set of states. The states with Ax = have large DOS 
near e = 0. For the state with A\ ^ these ZEBS are 
pushed to finite energies, away from e = 0. These shifts 
(splits) are largest for T> = 1 and decreases for small 
T> . As T> — + the DOS becomes independent of A\- 
We also calculated the junction energies for the different 
states. These are listed in Table GJ which shows explicitly 
that the Ax = 7r/2 states have lower energies. That the 
TRS-state cannot the minimum energy state and that 
the energy minimum state is at Ax = tt/2 (mod tt) may 
actually be expected from an argument based on symme- 
try and continuity. ||[l7) It is notable that in the small 
transmission limit, T> — 0.3, A\ = and A\ = tt/2 
have almost the same free energy. As seen in Fig. [j], 
the DOS for the two states are also similar except some 
small differences near e — 0. At this small T> the phase 
difference A\ ^ is inefficient in pushing the states that 
were orginally at e = to finite energies. 
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the structure at the interface. 
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FIG. 1. The magnitude and the phase \d of the order 
parameter r\d for a_t = and olr = 7r/4 at T — 0.2T C for 
the states corresponding to Ax = and to the energy min- 
imum with TRSB Ax / 0. The states with Ax / have 
a transverse current density j y along the boundary. In the 
lower panels c and d are the corresponding DOS on the two 
sides of the interface. Do is 1.0 in a and c and 0.3 in b and 
d. The units are fesTc for \r)d\ and 2e«/A i '/|r;(cxD)| for current 
densities in all graphs. £ = hvf /2nT c . 

The DOS recovers to its bulk value as one moves away 
from the interface. (Fig. ||) At x ~ 10£ o the bulk d-wave 
DOS is well-recovered showing only exponential tails of 
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FIG. 2. The spatial dependence of the DOS for an interface 
with T> — 0.7. The orientation is (q_l = 0, an = -|). The 
lower (upper) set of DOS are for the left (right) hand side of 
the interface. The DOS are sampled at a spacing of 1£ . The 
thick lines indicate the DOS at the interface location. 

The value of Ax where the interface energy is a min- 
imum depends on the orientations of the crystals, the 
transmission coefficient and temperature. An example 
is as shown in Fig | (c.f. Refs. @|||). The com- 
parison between the free energies for the orientation 
(aL,ctn) = (— 7r/12, 7r/6) is also shown in Table EL 
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FIG. 3. The phase difference Ax that minimizes the in- 
terface free energy as a function of oil- The misorientation 
angle 8 is kept fixed at 7r/4. The temperature is 0.2T C 

Recent experiments indicate that the oxide super- 
conductors probably also have an attractive s-wave chan- 
nel with a strength such that the bare T c for the s-wave is 
about 10% of that of the dominant d-wave. g] While in 
the bulk the order parameter is purely d wave, near the 
interface both d- and s-wave components can co-exist. In 
this case 



A(£/,ar) = ?jd(:r)\/2cos(2(<75 - a)) + r] s (x). 



(3) 



In Fig ^ we have plotted the order parameters r/^ , 77^ of 
the minimum energy states for (at, or) = (0,7r/4) and 
for different transparencies. For all D the phase differ- 
ence between the d-wave order parameters on the two 
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sides of the interface is 7r/2 as in the case without the 
s-wave component of the order parameter. In our gauge 
where the d-wave order parameter is real for x — > — oo, 
the s-wave component ry s is real for all x, being positive 
for x > and negative for x < 0. |Q The order parame- 
ter for x > is in the TRSB combination s + id. For the 
large transmission T> = 1 both the order parameter rjd 
and the DOS are qualitatively equal to the pure d x 2_ y 2 
state shown in Figs. |l|a and c. It is clear that the s-wave 
channel does not play an important role. It is rather the 
tails of the off-diagonal parts of the green's function of 
either side of the interface that are leaking into the oppo- 
site side that give the dominant TRSB. Hence, for large 
transmission the main mechanism for TRSB is the prox- 
imity effect even with a subdominant channel of mod- 
erate strength present. As V a is reduced the side with 
a = 7r/4 shows increasing pair-breaking due to reflection 
of quasiparticles by the interface and the TRSB gets more 
localized to this right hand side. This is also seen in the 
transverse current density, j y (x), which is much larger 
on this side. In the small T> limit the proximity effect is 
gradually shut off and the presence of the subdominant 
channel is largely responsible for the TRSB-state. 




-0.50 



-0.25 0.00 

e/2nT 



0.25 



±0.50 



-0.25 0.00 

e/2nT, 



0.25 0.50 



FIG. 4. The order parameters r/d and r/ s and the trans- 
verse current j y for oil = and ctR = 7r/4. T = 0.2T C and 
the sub-dominant T C 2 = 0.1T C . The transparency, T> , of the 
boundary is 1.0 and 0.3 in panels a and b. The corresponding 
DOS at the interfaces are in the lower panels c and d. 

The DOS shown in Figs. [j] and f| display considerable 
structure. These results are very different from those 
where the suppression of the order parameter near the 



interface is ignored (not shown) . In particular additional 
bound states are present at finite energies. These bound 
states are the result of Andreev-scattering processes due 
to amplitude changes in addition to sign changes in the 
order parameter. 

In conclusion we have investigated time reversal sym- 
metry breaking at interfaces, in particular those with 
high transmission. We have shown how this TRSB can 
be understood from the density of states and the free 
energy of the interface. 
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